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Abstract 

We present four-dimensional gauge theories in Minkowski spacetime which effectively 
generate in certain energy regimes five-dimensional warped geometries whereas, in gen- 
eral, the fifth dimension is latticized. After discussing in detail several general aspects 
in such theories we present a number of exactly solvable examples. We also point out 
how a particular case, defined in an A-sided polygon and having a Z^v symmetry, has a 
similar realization in an appropriate supersymmetric setting with D3-branes. 



1 Introduction 



A quite attractive idea that has received a lot of attention is that our spacetime is higher 
than four-dimensionaL This idea resonates well with string theory considerations, as 
well as with the more recent realization, in connection with a geometrical reformulation 
of the hierarchy problem, that such extra dimensions can be even macroscopic in size 
since this is allowed by the present day experiment evidence In most of the litera- 
ture extra dimensions is an assumption that one makes. Their number, size and other 
properties are constrained from theoretical consistency and experimental results. In such 
approaches the important point is that, certain phenomena in our four-dimensional world 
are just "projections" of Physics in higher dimensions, where the fundamental theories 
are defined. 

From a philosophical point of view it is quite different to demand that extra dimen- 
sions emerge dynamically and that the resulting theories are only effective descriptions of 
fundamental four- dimensional theories in appropriate energy regimes. Generically, such 
extra dimensions should emerge between certain energy scales, but the corresponding the- 
ories should look four- dimensional at very high energies in the ultraviolet (UV) and very 
low energies in the infrared (IR). That guarantees the four-dimensional nature of inter- 
actions at a macroscopic as well as at a microscopic level, provided the four- dimensional 
theories are well behaved, i.e. are renormalizable, asymptotically free etc. 

In studies of gauge theories at strong coupling using the AdS/CFT correspondence, 
we have constructed a model which shares the desired features we just mentioned 0. In 
particular, we constructed a D3-brane solution describing part of the Coulomb branch 
of TV = 4 SYM theory with gauge group SU (Nk) , in which vacuum expectation values 
(vev's) to two of the six scalars transforming in the adjoint of the gauge group were given. 
The scalar vev's were taken to lie in an A^-sided polygon and in this way the original 
5*0(6) global invariance of the theory breaks to an 5*0(4) x Z^. This theory turns out to 
look four- dimensional for very large and very low energies and effectively five- dimensional 
for a wide range of energies in between. This conclusion was also verified by examining 
the potential energy for a quark- ant iquark pair in [Q]. This potential, as we will also 
see in this paper, has the usual Coulombic bahaviour 1/r for very large and very low 
distances, but it behaves as in between, hence confirming the effective emergence of 
a fifth spatial dimension. 

More recently a similar idea was proposed in ||^ and in directly in a conventional 
field theory which need not be supersymmetric. A four-dimensional field theory was 
defined in an A^-sided polygon with gauge group SU{N) x SU{k). Fermions provide the 
necessary link and transform bilinearly under nearest neighbor pairs of gauge transfor- 
mations. If the dynamical energy scale for the SU{N) gauge theory A^v is much larger 
than that for the SU{k) theory A^, then for energies comparable to A^v the fermions 
condense in pairs since at these scales each of the SU (N) gauge groups becomes strong, 
whereas the gauge groups SU{k) can still be treated perturbatively. The authors of 
and P argued that the appropriate theoretical description is in terms of a latticized 
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gauge theory for SU{k) interacting with a a-model that provides the hnk between the 
different sites with nearest neighbor type of interaction. The gauge group is broken down 
to the diagonal subgroup SU{k) with the famihar Higgs mechanism. The vev's v were 
chosen to be the same at each site and are inversely proportional to the lattice spacing. 
It was subsequently shown that the theory looks four-dimensional for distances smaller 
than the lattice spacing 1/v and for distances larger than the size of the polygon ~ N/v. 
However, in between, for 1/f ^ r ^ N/ v the theory looks five-dimensional with the fifth 
dimension builded up by the Kaluza-Klein states of the large circle of radius N/v. The 
situation is the same as in ||^ in the sense that in the latter case the cr-model is provided 
by the six scalars in the A/" = 4 SYM theory that transform in the adjoint of the gauge 
group. Also in both cases the scalars assume vev's precisely the same way. This will be 
shown explicitly in section 4. 

We note that a similar idea of a latticized array of gauge theories has been proposed 
in the past in an attempt to explain the hierarchy of strengths between the electroweak 
and strong interactions 0. Also, there is a close relation to lattice gauge theories with 
anisotropic couplings in five or higher dimensions, which, for certain values of the cou- 
plings, exhibit four-dimensional layer-type phases as it was proposed in and further 
elaborated on, in @, |[. We also note recent related work that has appeared, with em- 
phasis into more phenomenological applications in |jTO|-|T^ and some connection with 
non-commutative geometry in [15, 1^. 

The purpose of this paper is to first provide in section 2, a systematic discussion of the 
general case where the latticized array is not a polygon and moreover the vev's depend on 
the site location (the case of varying gauge coupling in treated in the appendix). This will 
give rise, in certain energy regimes, to a five-dimensional description in a curved warped 
geometry instead of just the usual Kaluza-Klein case, where the fifth dimensional is flat. 
In section 3 we work out some examples completely in the lattice formulation, whereas 
for some others we resort to the continuum approximation. We believe that these results 
can be further used towards constructions that are more phenomenologically oriented. 
Finally, in section 4, we review and refine certain aspects of the model introduced in in 
connection with the more recent developments. We also suggest a natural mechanism for 
the resolution of singularities in solutions that appear in studies of the Coulomb branch 
of the Af = 4 SYM theory. We have also written an appendix where we treat the case 
of varying gauge coupling according to the lattice position (but constant vev). We show 
that this is related to the case of varying vev's and constant gauge coupling of section 2, 
by a generalization of the formulation of the usual supersymmetric quantum mechanics, 
but now on the discrete lattice. 
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2 The general formalism 



As in Q we start with an effective four-dimensional low energy action which is the sum 
of YM actions for SU{k) as well as for a Higgs field that provides the necessary linking 



^ „ TV N-l 

^ = ~ 4(iV+l)^| J d'^Y.J'(^P' + Y.tT{{D,^,yD^^,) , /x = 0,l,2,3. (1) 

Each gauge term is labeled by an integer j, whereas the Higgs field $j links the different 
gauge terms. In transforms in the (Nj,Nj+i) representation of the gauge group, where 
we note that group generators corresponding to different sites commute. Accordingly the 
covariant derivative is defined as D^^^j = d^^j — iA^^^j + i^jAjj^^. The gauge coupling 
{N+l)gl is associated with each factor SU {k) separately. On the other hand the coupling 
gl is associated with the diagonal SU{k) that is left unbroken and dictates the physics at 
very low energies. We assume that there is a potential for the Higgs field $j that allows 
this field to develop a vev Vj = vfj, where we find it convenient to separate an energy 
scale V from a set of numbers fj that characterize the strength of the vev at each site j. 
This is our main difference from 0, Q , which however will have important consequences. 
A similar action to ([l|) has appeared before in section 2 of where in a five-dimensional 
f/(l) gauge theory on an anisotropic lattice, the continuum limit in the four- dimensional 
layer was taken, but the fifth dimension was left latticized. 

We digress, to note that, the case with varying gauge coupling, and constant vev's 
are completely related to the case of constant gauge coupling and varying vev's that we 
examine here. This will be shown in the appendix. 

Turning back to (0) we see that the relevant part of the action corresponding to the 
gauge field fluctuations is 

^ „ N N-l 

The mass term above can be written in terms of an (A^-|-l) x (A^-|-l) matrix T defined 
by writing it in the form A^^TA^. Finding the mass spectrum amounts to simply diago- 
nalizing the matrix T. However, this can be a cumbersome procedure except for a small 
number of sites. In order to develop a systematic approach we prefer to write this mass 
term as 

N 
j=0 

where the operator Tj is defined as 

T,. = _ ^)^2(gd, _ ^) ^ _4^2 sinh(rf,/2)/^_ , sinh(rf,/2) , (4) 

with dj being the usual differential operator satisfying [dj,j] = 1. The rewriting of 
the mass term is correct up to a term which is taken care of by the boundary conditions 
as discussed below. 
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In principle, in order to find the spectrum of tfie gauge field fluctuations we need to 
diagonalize the operator Tj. This can be done in some cases explicitly, but of course not 
in general. Before we come to specific examples it is important to mention some generic 
properties and features. Consider the eigenvalue equation 

Tjcpj = M^^j . (5) 

By construction the operator Tj is Hermitian and therefore it has real eigenvalues. In 
addition, since it can be written in the form Tj = qjqj, where qj = vfj{e'^^ — 1), its 
spectrum is positive semi-definite , i.e. > 0. The eigenvalue equation (^) should be 
supplemented with some boundary condition. In order to find what consistent boundary 
condition can be imposed we first realize that the equation of motion for as it follows 
with the mass term (^) is different than that following with a mass term as in (^, by the 
extra term v%{A'^'^^ — A^). Translating this into a condition for the eigenvectors (pj we 
have 

boundary condition I : (pN+i = 4>n ■ (6) 

This boundary condition preserves the zero mode. An alternative boundary condition 
which does not preserve the zero mode is 

boundary condition II : 0o = and (or) (j)^ = . (7) 

The physical motivation and consistency for such a boundary condition needs some fur- 
ther explanation. Let us relax the assumption that the gauge coupling is the same at all 
sites and introduce a "defect" by assuming that one of gauge couplings at the site j = N 
is much smaller than the couplings at the other sites which are taken to be equal to gf, 
i.e. gl j^ <^ g\. Hence, in this limit, after we scale A^ — >• A^g^^^, we see that the gauge 
field A^ drops out of the mass term in (Q). In terms of the corresponding eigenfunction 
0Ar this implies the condition = in (j^). A similar comment applies if we take the 
gauge coupling at the site j = very small, which implies then the boundary condition 

00 = 0. 

Imposing the boundary condition we find that there is one massless eigenvalue 
corresponding to the zero model, as well as massive ones. We emphasize that this 
procedure of finding the eigenvalues is completely equivalent to simply diagonalizing the 
(A^ + 1) X (A^+1) matrix T defined before by writing the mass term in (Q) in the form 
A^^TA^. We can always construct a set of A^+1 eigenvectors {0j>, n = 0,1, . . . , N} of 
the operator Tj that obey the eigenvalue equation (|]) as well as the orthogonality and 
completeness relations 

N N 

Yl ^h^^lrn = ^n,m , Yl ^h^^l,n = ^j,k ■ (8) 

j=0 n=0 

The normalized wavefunction corresponding to the zero mode is (pj^ = {N + 1)"^/^ with 
Mn = 0. 
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On the other hand, imposing one of the boundary condition in we find that there 
are N massive eigenvalues, but the zero mode is projected out of the spectrum. Then, 
in the finite sums in (|^) we exclude the corresponding terms. 

Lets us introduce a source term corresponding to a unit charge at the origin of the 
coordinate system and at the k-th site. Then, the differential equation that determines 
the potential energy at site j due to a charge located at site k is given by 

VlVj^k - Wj^u = 47r(iV+l)(725^.,<5(3)(x) . (9) 

We may expand Vj is terms of the basis eigenvectors of the differential operator Tj as 
^j,k{r) = Y^n Vn,k{^)4>j,n- Usiug (D we eventually arrive at 

= -^^^^ f:0„.0^,.e-^- . (10) 

n=0 

We may discuss quite generally]^ several features for the behaviour of the potential 
in ([T0|). For instance, at distances much larger that the typical size of the array, i.e. 
r ^ N/v, obviously only the zero mode, if it exists, contributes to the sum in (p!0|). 
Hence we have that 

Zero mode exists : ^■,fc('") — — - , for ^ ^ N/v . (13) 

Therefore, at low energies we obtain the usual Coulombic behaviour with coupling con- 
stant gl (we ignore possible slow variations due to the running of couplings). This is the 
expected behaviour for a four-dimensional gauge theory for SU{k) which is the unbroken 
diagonal subgroup of the original 1 copies of SU{k). This behaviour is universal for 
all sites due to the fact that the wavefunction for the zero mode has the same value at 
all sites. For the cases that the zero mode does not exist the typical behaviour is 

Zero mode does not exist : Vj Jr) ~ _ ^40j,i0A:,i ^-A/iT- ^ iov r ^ N/v , (14) 

r 

where the inverse of the lowest eigenvalue, which is expected to be comparable to the 
typical size of the array, i.e. Mf^ ~ N/v, sets the range of the Yukawa interaction. Its 

^For completeness, we note that, when the eigenvectors (fij^n form a complete but not orthonormal 
set, we may follow the procedure of Gram-Schmidt in order to construct such a set. Otherwise, (H) is 
replaced by 

N N 
j—0 n,m— 

where Amn is some (A^+1) x (A^+1) Hermitian matrix. Consequently, (|lo| ) is replaced by 

V,Ar) = -^^^^±1M f2 0,,«^„^0^.„e~^^- . (12) 

n,m— 
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effective overall strength depends on the sites.0 In contrast, for distances r ^ 1/v all 
modes in the sum in (^) contribute equally. After using the completeness relation in (^) 
we obtain 

^,>(r)^- ^^^^^^^' (5,,fc, as r^l/v. (15) 

This is the expected behaviour for a four-dimensional theory at high energies where we 
have + 1 independent copies of gauge theories for SU{k), each with gauge coupling 
{N+l)gl. Since, the various copies are independent the potential is non- vanishing only 
for charges charged with the same SU{k) factor, i.e. j = k. For charges charged with 
respect to different SU {k) factors where j ^ k, the potential has an expansion in powers 
of r with leading term of order g\vN for site-separation of order 1, and g\v/N for site- 
separation of order N . Hence, the magnitude of the force between charges located at 
the same site is much larger than that between charges located at different sites. In the 
latter cases, higher loop diagrams that contribute may be considered as, depending on 
the value of the coupling constant, they may give comparable contributions (a similar 
comment has been made for a specific example in |^). 

We conclude our general discussion with further remarks on the general behaviour 
of the potential as a function of r. First, notice that, as it follows from (0), we have 
implicitly chosen opposite charges for the two charged particles with respect to the diag- 
onal SU [k) subgroup that prevails at low energies. The question is whether it is possible 
in some models to have a repulsive instead of attractive force for some ~ 1, although 
of course eventually the force has to turn attractive according to (P^. Certainly, the 
possibility of an repulsive force for some range of r is excluded for charges at the same 
site since all terms in ([To|) come with the same sign when j = k. However, for j ^ k 
not all terms have the same sign and the possibility of a repulsive force is not excluded 
for some range of the values for r. Nevertheless, in all of our examples in this paper the 
force remains attractive for all r, in agreement with physical expectations. 



2.1 Continuum limit and a fifth dimension 

A quite interesting behaviour occurs for intermediate distances, which are much larger 
that the lattice spacing, but also much smaller than its typical size ~ N/v, namely for 
1/f ^ r ^ N/v. Here we assume that A^ ^ 1. Then we may approximate the lattice 
by a continuum with fifth space- variable = j/v. Then the differential operator Tj is 
well approximated by T^x^) = —d^f{x^Yd^, where f{x^) is the continuum limit of the 
sequence of numbers {fj} which is assumed to exist. Then, the continuous limit of the 
action is written as 

^ = / rf'^^^5tr {Fl + f{x,)\d,A^f) , gl = Ngl/v . (16) 

^In ( |r^ the strength of the Yukawa interaction is constant. This is due to the discreteness of the 
spectrum. In cases where the continuum hmit is considered this strength is r-dependent, giving rise to 
an overaU factor for e"*^^*" which is proportional to instead of 1/r. 
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This can be reinterpreted as the five- dimensional action —-^ J d!' x\/ —Gti{Flj ^) , where 
M = {fi, 5) and in the gauge = 0, defined in the background with a metric given by 
the fine element 

ds'^ = fix^frif.^dx^dx" + dxl . (17) 

Hence, what we have derived is that there is an energy regime where we may safely 
use a five-dimensional gauge theory defined in a curved, in general, warped spacetime, 
in order to compute physical processes in a gauge theory in the usual four- dimensional 
Minkowski spacetime. The four-dimensional gauge theory is renormalizable whereas the 
five-dimensional one is not and only emerges as an effective description in the appropriate 
energy regime. 

We make the important note that the regime where 1 ^ f r ^ is a strict continuum 
limit where the operator Tj is well approximated by the second order differential operator 
T(x5) = —d^f{x5)'^d^. Higher terms that appear in the derivative expansion of Tj are 
neglected. However, there is a less restrictive continuum limit that retains all such higher 
derivative terms. As it will also be confirmed in the examples below, this limit is attained 
for vr <^ N and therefore is valid for distances tr ~ 1, where the approximation that 
leads to (|T6|) fails. 

For future convenience it is better to use the variable z instead of x^, with change 
of variables as dx^ = e^dz. Also we will rename by e^^. The equation determining 
the potential between two charges is still given by (P) with Tj replaced by T^x^). The 
corresponding eigenvalue equation for this operator becomes 

T0„ = -e-^d,e^dMz) = Mlct>n{z) , (18) 

where we have used the variable z instead of x^. The orthogonality and completeness 
relations read 

I dze^M^Wn^i^) = 8n,m , J] M^Wni^') = C'^Siz-z') . (19) 

n=0 

As in the discrete case, we may distinguish the following two cases: If the spectrum of 
the operator T has a mass gap then the behaviour of the potential between charges is, 
for large distances, of the Yukawa-type, whereas if there is no such mass-gap we have a 
power-law type of behaviour. 

We note that, even though for notational purposes we have used the discrete set of 
number {n G Z} in order to label our states, the spectrum can be discrete of continuous 
depending on the function A{z). In order to check that it is convenient to cast the 
differential equation in (|^) as a Schrodinger differential equation. Indeed, after defining 
(pn = e~^/^\E'„ we obtain 

- dl^n + ^Sch.^n = M^^, , l^sch. = ^(5.^)' + ^O^A . (20) 

This potential has the form of the potentials appearing in supersymmetric quantum 
mechanics Jl^ with corresponding superpotential W = —^dzA. Hence, the spectrum 
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is positive semi-definite as already expected from the fact that it originates from the 
spectrum of the operator Tj which we have already argued that it is positive semi- definite. 

If the Schrondinger potential in (^) can be neglected for some wide range of values 
of the variable z, then in this range the wavefunctions are approximated by plane- waves 
as '^k — l/V^e'^''^- Then, using that (pk = e~^/'^'^k = /(a^s)"^^^^/; we obtain that the 
potential between charges at positions z and z' in the fifth dimension, is given by 

2 ^^(A(z)+A{z') 

y^Ar)^-^--Y-7 ^- (21) 

TT + [z — z')^ 

This behaviour resembles the usual Kaluza-Klein case, but with an effective coupling 
that depends on the fifth dimension. This is due to the fact that the effective fifth 
dimension that emerges is not fiat. Such behaviour will be encountered in some of the 
explicit examples below. 



3 Solvable examples 

The examples that we specifically work out in this section belong to two categories: 
Those where we work directly in the lattice formulation and those where the continuous 
approximation based on (0) is employed from the very beginning. We first consider the 
case of a periodic array on a polygon that has been considered before in this context in 
IQ, ^ and its supersymmetric variant already within the context of the AdS/CFT cor- 
respondence in 1^. In our second example the periodicity condition is relaxed, resulting 
into a model that has a discrete mass spectrum, but no zero mode. Both cases share the 
feature that the vev's are the same at all sites. Our third example has vev's that vary 
with the site position, but even so it turns out that is is completely solvable with the 
help of Laguerre polynomials. Finally, within the continuous approximation we consider 
two classes of examples with varying vev's. Both have a continuous spectrum but one of 
them has a mass gap. 



3.1 Periodic array on a polygon 

In this example we give the same vev at all sites, i.e. vj = v. Moreover, because of the 
Ztv+i symmetry of the array, we require periodicity of the wave function as 0j+Ar+i = 
Since we have that fj = 1 the operator in @j is simply Tj = —Av'^smh.'^{dj/2). 
This commutes with the differential operator dj and therefore they have the same set of 
eigenfunctions which when properly normalized read 

= ^=^6^-^ , j = 0,l,...,iV, n = 0,l,...,iV. (22) 
The N + 1 mass eigenvalues are given by 

M^ = 2vsmlj^] , n = 0,l,...,N, (23) 
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where we also observe that a zero mode exists. Since |0j>| = \/\/N+\ the potential (|TOD 
is the same for all sites, as expected due to the Z^j^i symmetry of the array. Even in this 
simple case, it does not seem possible to perform the sum in ([T0|) explicitly.^ Nevertheless, 
the behaviour for r ^ N/v and r -C 1/f is that predicted on general grounds by (|13|) 
and (|T5|), respectively. 

The continuum limit holds for vr ^ N and j — k ^ A^, Then the original discrete 
periodic array is replaced by a continuum periodic array, but we retain all terms in the 
derivative expansion of the operator Tj. The expression for the potential turns out to be 
given by the following integral 

V,,k{r) = -2^ / dv;cos(2(j-A;)^)e-2--'^^ . (24) 
Jo 

We may evaluate this in terms of generalized hypergeometric functions, but we will refrain 
from presenting the result since the resulting expressions are quite complicated. 

If we further restrict to distances 1 -C f r -C we find the behaviour 

V,,,(r) ^ ^ , for l/v«r«N/v, (25) 

where 0:5 = j/v and x'^ = k/v. This is the regime where the description in terms of 
a five- dimensional theory in flat space (since f{x^) = 1) becomes appropriate. Then 
the continuum periodic array is replaced by a an infinite continuum array. In that case 
the properly normalized wavefunctions are (pki^^) = e^'^^-' / \/2it, with z taking values in 
the whole real line, whereas the mass spectrum = fc^ is continuous. Then using that 
V = —gl/r dke'^'^^'^ / {2n) we exactly reproduce (pSj). This is the typical Kaluza-Klein 
behaviour corresponding to the five- dimensional effective theory (p!6D and was already 
found in [g g. 

There is in addition the scale regime vr ^ 1 where the original discrete periodic array 
is replaced by a infinite, but nevertheless discrete array. Then the potential is given by 
an infinite sum which can be computed exactly 

2 00 00 ^ 

94 \ ^ 2ni{j-k)n/N -2nvr\n\/N _ _ 2 



V-,(r) - -ii 27riU-k)n/N -2nvr\n\/N _ _ 2 ST^ 

''"^ ' - r ^ - 2^ ^2 + ( _^_jVn)7t;2 

n=— oo n=— oo ' 

gj sinh(27rrf/A^) 

(26) 



r cosh(2'7rrf/A^) — cos (27r(j — A;)/A^) 
This expression interpolates between those in (p5D for 1 ^ f r -C and ([131) for vr ^ A^. 



In the former case the corrections to (^Sf ) are in powers of vr/N, whereas in the latter 
case the corrections to (0) are Yukawa-like, corresponding to the Kaluza-Klein modes 
of a reduction to a circle of radius N/v. 

^However, we may easily prove that for r ^ l/w, the behaviour of the effective potential is Vj^k — 

-6j,k94{N+l)/r-TTvgl/{N + l) (sin^ {TT{j - k)/{N+l)) - (i7r/(A^+l))^) V---. This is in accordance 
with the general behaviour that we have stated after eq. ([l5|). 
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3.2 Relaxing the periodicity condition 

As in the previous example we give the same vev to all sites, but we no longer impose 
the periodicity condition. Instead, we impose the boundary condition (|^) that the wave 
function (f)j vanish at the endpoints, namely that 0o = 0Af = 0. The properly normalized 
wavefunctions that form a complete set are 



— sm{7rjn/N) 



J = 0,1,...,N 



n = l,2,...,iV-l . 



The mass spectrum is given by 



/ Tcn\ 

Mn = 2t; sin — 



n = 1,2 N-1 



(27) 



(2^ 



where the zero mode is projected out. With these eigenfunctions and eigenvalues the sum 
in ( [TDI ) cannot be computed exactly. Nevertheless, we may show that the behaviour for 
r ^ N/v and r ^ l/t> is that predicted on general grounds by ( [TB| ) and (|TB|), respectively. 

In the limit vr ^ 1 and for ^ 1 we have that (we set L = N/v) 
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sin(7rjr;,/A^) sin(7rA;n/A^)e' 



■■Kvnr/N 



n=l 



Sin (f ) sin (f ) sinh (f ) 



cosh ijjf) — cos 



N 



cosh (^) — cos 



N 



(29) 



For r <^ X5 ~ Z5 <^ L, where = j/v and X5 = k/v, we obtain the same behaviour as 
in dl). Also 

^^'^"^^"^^ , X5, x's « r « L , (30) 



V ' ~ 

•^5 5*^5 



which is a behaviour similar to that of a seven- dimensional theory. Finally, for distances 
much larger than the length L we have 



2gl 



sm{7cj/N) sin(7rA;/A^)e-^^/^ , for r > L 



(31) 



in accordance with the general formulae ([TB|). 

There is also a continuum limit for distances rv -C A^, that goes beyond the two 
derivative approximation (|l^ 



vrr 



7r/2 



d(psm{2j{p) siia{2k(p) e 



-2vr sin ip 



(32) 



Again this is expressible in terms of generalized hyper geometric functions. 
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3.3 Array with varying vev's 



Cases with vev's varying according to the site position are in general more difficult to 
handle for increasing site number A^. However, let us consider the case of Vj = f V7; 
j > 1.Q Though not obvious, it is possible to completely determine in general the 
eigenvalues and eigenfunctions of the operator Tj which now equals 

Tj = v\e-'^^ - l)j(e'^^ - 1) . (33) 

As in some quantum mechanical problems it is convenient to pass to the momentum 
representation where the operator dj is just a number denoted by ip and the site posi- 
tion j becomes the differential operator idp. Then the eigenfunctions in the momentum 
representation satisfy the differential equation 

tv^e'P - l)dp {{e'P - l)0p) = , (34) 

with properly normalized solution^ 



t e 



im^ /2cot{p/2) 



e*'"'/2=°*(p/2)(l+icot(p/2)) 



^ 0r e'P-1 2^1 

< p < 27r , m = Mjv . (35) 

In order to compute the potential between charges we need the wavefunction in the 
position representation This is given by 

^ J dpsin{jp)(f)p . (36) 



After changing variables as a; = cot(j9/2) we have to compute the integral in (|36D which 
takes the form 

The integrand has poles at x = ±i in the complex plane and we may compute it using 
standard methods in complex analysis. Since > we may choose the contour of 
integration to close in the upper half plane and to include x = i. Then, only the first 
term in the integrand in (|37|) contributes. The result is 

^ {] - 1)! (9xJ-i L ^ ' 

^(j-l-n)!(n!)2 • ^ 



"'For convenience, we have shifted the range of j by one. Hence, j starts from j — 1 and not fr'om 
j = as in the general formulation before. We take this into account in the rest of this subsection. 

^We note that the case of an array with Vj = vj, j > 1 is also solvable by applying the same 
idea. Indeed, in this case after we substitute 0p = (1 — q)'-'-F{q), where q — cos^(p/2) and /i = 1/4 + 
j/2yTn^ — 1/4, we obtain that F(q) satisfies a hypergeometric equation with a — b — ^ and c — 1/2 in 
the standard notation. We also note that, normalizability of (f>p requires that the mass spectrum has a 
gap, Mgap = v/2. 
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One recognizes the last sum as an expression for tlie Laguerre polynomials Lj_i(m^). 
Therefore 

0. = e— K) , j = l,2,... . (39) 

Having obtained this solution in a systematic way, it is important to verify directly that 
(^) indeed satisfies the eigenvalue equation Tjcpj = M'^(j)j, where Tj is given by (|33|) . 
This is quite straightforward, using the recursion relations for the Laguerre polynomials. 



3.3.1 The limit of infinite array 

In the case that the array is large, so that ^ 1, boundary effects can be neglected for 
vr <^ N. Then the mass spectrum is continuous with no mass gap (M^ > 0). Using the 
orthogonality and completeness properties of the Laguerre poljTiomials for G (0, oo) 
we see that (||) are satisfied. The potential between charges at sites j and k is then 

V,-fc(r) = -2^/ c/mme-'"'-'"^"Lj_i(m2)Lfc_i(m2) , for vr <^ N . (40) 
''^ Jo 

For specific values for j and k we may give the result in terms of the error function 
erf(fr/2) and powers of vr, but we will not present any results explicitly. 

It is easy to see that for f r <^ 1 we obtain (|I5[) whereas for vr ^ 1 and small values 
for j and k we obtain an fall-off for the potential as 

V.kir) ~ -2^ , for j-k<^vr <^N . (41) 

This is characteristic of the behaviour of a six-dimensional theory. 

The limit of site position j ^ 1 and ^ 1 is also of particular interest. Then using 
properties of the Laguerre polynomials we find that 

Jo (2m a/7) , as j ^ 1 , ^ 1 and jm? = finite , (42) 



where Jn is the zeroth order Bessel function. This is the result that one would have 



obtained from the five-dimensional two-derivative effective theory based on (0), as we 
will shortly verify (see footnote 6 below). 



3.3.2 Finite arrays 

In the case that the array is finite with j = 1, 2, . . . , + 1 we have to consider boundary 
effects. There are two choices, one projecting out the zero mode and one preserving it. 

Projecting out the zero mode: We impose the boundary condition (|^) that the wavefunc- 
tion vanishes at the end, namely that (pN+i = 0. This condition determines the allowed 
values for as the zeros of the corresponding Laguerre polynomial Ljv 

L^(m2) = 0, n = l,2,...,iV . (43) 
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The resulting eigenfunctions (j)j^n = e~'""/^Lj_i(m^), j = 1,2, ...,A^ with n = 
1,2, . . . , N form a complete orthogonal set, but there are no longer properly normal- 
ized to 1. Using the Darboux-Christoffel formula for the Laguerre polynomials we find 
that the correctly normalized orthogonal complete set of eigenvectors is given by 



L,.,{ml), 3 = 1,2,. ..,N, n = l,2,...,N. (44) 



Then the potential between charges located at the lattice-sites j and k, and separated 
by a distance r is 

V^r) = -^i^^ L,_iK)L._,(mDe--- . (45) 

As an example consider the simplest non-trivial case with N = 2. The two roots of 
L2(m^) = are given by = 2 =f ^/2. Using the fact that Lq{x) = 1 and Li{x) = 1 — x 
we find that the two eigenfunctions forming a complete orthonormal set according to (|^), 
are 

The potential between charges contains of course only Yukawa-type terms. We have 

V^,ir) = -M(a-.e-" + aje-+-) , ^,J = l,2, (47) 
where the 2x2 matrices are given by 

a^ = ?f"f ^^'). (48) 



4 V ±V2 2tV2 

Keeping the zero mode: Let us now impose the boundary condition (||), that is (pN+2 = 
4>N+i, or using the properties of the Laguerre polynomials, 0'^+2 — 0- That determines 
nonzero mass eigenvalues m^, n = 1,2, ... ,N from the condition 

L'N+ii^l)=0 ^ LN+iiml) = LNiml) = L'^iml) , n = l,2,...,N. (49) 

However, there is also the zero mode which is preserved by the boundary condition. As 
before, using the Darboux-Christoffel formula for the Laguerre polynomials we find that 
the correctly normalized orthogonal complete set of eigenvectors is given by 

1 



Lj-i{ml) , j = l,2,...,N+l , n=l,2,...,N 



VN+l\LN{m^ 
1 



J,n I , , „^| ^] — LV"'nl 1 J x,^,...,-it I J- ; j.,A^,...,±y , 



3 = 1,2,. ..,N + l. (50) 



Then the potential between charges at sites j and k has a Coulombic as well as Yukawa 
terms 

2 2 ^ 1 

= -7 - 7 E h_Ml)Lu^r{<)e--- • (51) 
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As an example consider again the case of = 2. Now there are three mass eigenvalues 
given by m-Q = and the two zeroes of the algebraic equation L'^{m'^) = which are 
= 3 =F ^/3. The corresponding eigenf unctions are 

The potential between charges is 



2 2 

= ~ ^^^Jk^~" + ^%^"^^'"') ' J, = 1, 2, 3 , (53) 
where the 3x3 matrices are given by 

2±v^ -1 -It V3 
b^ = - -1 2^73 ±v^-l I . (54) 



2 



-l^v^ ±V3-1 2 



Finally, we note that when becomes very large, the spectrum in both cases that 
we have considered becomes continuous and the above expressions approach the corre- 
sponding ones in ( |39D and (^). This can be seen using the property of the Laguerre 
polynomials 

LJx) = —=-. — --rr cos(2^/nx - IT / 4) + 0(n~^^^) , for n > 1 . (55) 



3.4 Example within the continuous approximation 

We consider now cases where the lattice is already treated in the continuum approx- 
imation where a description in terms of a gauge theory in a five-dimensional warped 
spacetime (|16|) is valid, namely when 1 -C tr <^ A^. The distribution of vev's might be 
such that the mass parameter that sets the scale at these distances is not v itself, but 
some other mass parameter Mq. This is assumed to be such that v/N <^ Mq <^ v which 
implies that the corresponding length scale I/Mq is within the allowed distances for the 
continuum approximation to be valid. 

First, let us consider the case with /(xs) = 2M0X5, xt^ > 0. Applying the general 
formalism we easily find that the variable 0:5 and the variable z that enters in (|18]) are 
related by 2MqX5 = e^*^°^, with —00 < z < 00 and that the function A{z) = 2Mqz. 
Therefore the Schrodinger potential in (pO]) is just a constant, i.e. Vsch. = Mq. That 
defines a mass gap since in order to have normalizable wavefunctions we must require 
that > M^. We obtain 

„ — M()2 

Mz) = ^{M' - Mo2)-V4e^(*^^-*^o)^/^^ . (56) 

2a/ TT 

Hence, the potential between charges at positions z and z' in the fifth dimension is 

V.,.' = - ' ° ^ K.iMoR), R=sjr^ + {z-z')\ (57) 
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where Ki is the modified Bessel function of order one. If we consider distances much 
smaller than the characteristic scale of the mass gap I/Mq, then we obtain indeed the 
behaviour (pT]). Instead, in the opposite limit of distances much larger than 1/Mq we 
find that 

K,.'(r)^-^e-^^»(^+^')— -^e-^^«^, for 1/Mo « i? « iV/^ . (58) 

Naturally, the range of Yukawa interaction is set by the mass gap, whereas the strength 
has the Xjr^l'^ behaviour as noted in footnote 2. 

Also, the value of the mass gap Mq is in agreement with what was noted in footnote 
5 for the exact lattice treatment of the vev distribution Vj = vj, after we set Mq = v/2. 



3.5 Another example within the continuous approximation 

A generalization of the previous case which however does not lead to a mass gap for the 
spectrum, is to take /{x^) = (2(z/ + l)MoX5) ^("+1) with 0:5 > and where the exponent 
of X5 is introduced as such for later notational convenience. We find that the variables X5 
and z are related as 2(z/ + 1)1/0X5 = {Moz)^^''+^^ and also that = (Moz)2^+^ Hence, 
the Scrhodinger potential takes the form 

Vsct. = ^ , 59 

from which we conclude that the mass spectrum is continuous without a mass gap. The 
eigenfunctions are given in terms of Bessel functions as 

0m(^) = -^^(Mozy'^MMz) . (60) 

It turns out that^ the potential between charges at the same position z in the fifth 
dimension is computed in terms of a hypergeometric function as (for simplicity we do 
not deal with the case of charges at different positions in the fifth dimension) 

^ -?mm^Jl^F Urn..- V, 2. . 1 



0Fr(z/ + l)r2(r2 + 4z2)^'+i/2 V ' + Az^ ^ 

for l/i; < r -C N/v . (61) 

We also note the limiting case 

and that for r <^ 2; we confirm the behaviour (PT|). 

^Note that, for Mq = v/2 and = we have that Vj = w^/J. Then indeed the wavcfunction ( |60| ) 
agrees with (E2h. 
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4 Strong coupling considerations 



We would like to study similar effects at strong coupling using the AdS/CFT correspon- 



dence |]I8| . We will see that our findings and the emergence of an effective fifth dimension 
are in precise agreement with the gauge theory considerations so far. 

The clearest discussion for our purposes can be done for the Coulomb branch of 
the Af = 4 SYM theory with gauge group SU{Nk). This theory has six scalars , 
/ = 1, 2, . . . , 6 in the adjoint of the gauge group. We will reanalyze and refine the model 
introduced in and we will study the potential between a quark and an antiquark. 
We give vev's to two of the six scalars in the theory in a way that the original 5*0 (6) 
global invariance of the theory that rotates the six scalars, breaks to 5*0(4) x Zi<^. We 
may represent these scalars as Nk x Nk traceless matrices. In general, since the scalar 
potential in the SYM action is proportional to tr[X'^, X"^]^ we can choose the vev's 

along the diagonal so that the potential is zero. That preserves sixteen supercharges and 
breaks the conformal invariance. We arrange the Nk x Nk matrices X^'s in the form 
of block-diagonal matrices with N diagonal k x k blocks, where in each block the vev's 
are the same. The J\f = 4 SYM theory action contains the term tr(Z}^$^)^ which 
couples the gauge fields and the six scalars. It is convenient to choose the standard real 
basis for the SU{Nk) generators Jy = eij — -^Sijl^kxNk^ where the matrix elements of 
the matrices e^j are {eij)mn = ^im^jn- Hence, we give vev's to the scalars represented 
by the six-dimensional vector $, as $ = hiX^^^, where hi = J a are the diagonal Cartan 
generators. The masses of the gauge fields arise from the term £mass ~ X]/=i tr[$'^, A^'^ . 
After we define = A^^Jij we have 

N-l 

Cm^ss ~ k J2 (Ke.- XL - XL- KJ^A^'' ■ (63) 

i,j=0 

This basic formula can now be applied to specific cases. As in we choose for the 
vev's in the jth block the values 

= (0, 0, 0, 0, ro cos(27rj7Ar), ro sin(27rj7iV)) , j = 0, 1, . . . , X-l . (64) 

These form an X-sided polygon enclosed by a circle in the x^-xq plane of radius vq. Then 
we find that (0) becomes 

N-l 

Cm^s ~ r',k sin'(7r(j - z)/N)A'^A>^'^' . (65) 

i,j=0 

From this we read off the masses as 

M„ ~ roV^sin(7™/X) , n = 0,1, . . . , N -1 , (66) 

which are the same as the ones we found in (pSf). They have a degeneracy dn which is: 
do = N-l and dn = 2{N-n), for n > 1 (Note that 'EnZl dn = iV^-l). Using the 
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fact that Mn = Mi^_n it is easily seen that the expression for the potential is given by 
V = -gl/{Nr)Y^^Zldne-^'^"'- = -gj/r (j2nZle-^^"' - l/^)> which, for large enough 
A^, is the same as the potential ( p!OD corresponding to the mass formulae in (p3D (for test 
particles at the same lattice cite). That makes apparent that the model in p[ and the 
construction of (reviewed in subsection 3.1) describe essentially the same Physics. 

The advantage of having a D3-brane realization in a model with an emerging latticed 
fifth dimension, is that we may consider the theory at strong coupling using the AdS/CFT 
correspondence. In the type-IIB theory we consider a large number Nk of D3-branes that 
are distributed in the transverse to the branes six-dimensional space according to (p^). 
Namely, the gravitational metric that we will consider is 

ds'^ = H-^/^r^^^dy^dy" + H^^^dxl + ■■■ + dxl) , (67) 

where y'^, /x = 0, 1,2,3 are the parallel to the brane directions. The harmonic function 
H in the six-dimensional space spanned by Xj, i = 1, 2, . . . , 6 is given by 

N-l 4 

H = y (68) 

where Rf = Anglk is the 't Hooft coupling associated with the SU {k) gauge theory. We 
see that the vev's in the field theory ( |64D become the centers in the harmonic function 
above. We emphasize that the supergravity approximation is valid even close to each one 
of the centers if we choose that 1 and also small values for g\ which implies that 

there are many branes located at each center, i.e. A; ^ 1. The number of such centers 
need not be large. 

The harmonic function H has been computed in [Q. It takes the form 

H = ""-—^ — E;v , (69) 

where i?^ = AnglkN. Also U"^ = xf + . . . + xj, u"^ = xl + and 

^ sinhA^x ^ ^(([/^ + Tq)^ — 4ro'U^)^^^ cosh A^x cos iV-?/; — 1 

cosh A^x — cos Nip + (cosh A^x — cos Nipy ' 

where ip is the polar angle in the xs-xe plane and 

e..i^+,/f^E5yr;. (n, 

Notice that H explicitly processes the symmetry 5*0(4) in the four- dimensional subspace 
spanned by Xi, X2, X3 and X4 as well as the symmetry under the shift ip — >• ip + 27i/N. 

The potential between a pair of heavy quark and antiquark separated by a distance 
L is foundQ by minimizing the Nambu-Goto action of a fundamental string in the su- 
pergravity background (^). We will not give details as the necessary techniques we 

^We denote the distance between the charges by L in order to confirm with standard notation in the 
AdS/CFT Uterature. This is nothing but the distance denoted by r in the previous part of the paper. 
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first developed for the conformal case in |T^. For backgrounds of the type that we have 
here we will use eqs. (5)- (7) in ||^. The trajectory that minimizes the action lies on 
the x^-Xq plane where u = U . However, not every angle if) is allowed, but only those 
with sm.{Nip) = 0. All angles satisfying this condition are equivalent because of the Z^- 
symmetry. For concreteness we may choose the trajectory along the xs-axis (for ip = 
and u = U). Let's define the function / = R^/H. Then for the trajectory that we are 
interested in we have 



fiU) 
E 



(f//ro)^ - 



S 
1 



+ 



2N 



(72) 



where we have used the fact that = t//ro, when u = U and ip = Q. Then the integrals 
for the length L and the energy V are given by 



f{u)U{u)-f{u,)) ' 



and 



V 



IT 



dU 



Uo 



fiU) 



f{U) - f{Uo 



1 



ro) 



(73) 



(74) 



where Uq denotes the lowest value of U that can be reached by the string geodesic. In 
principle, we have to evaluate the integrals and try to express V" as a function of L. This 
is impossible in general though it is easy to deduce that is a monotonously increasing 
functions of L from —oo to in the interval L G (0, oo). 



4.1 Examining various energy regimes 

As usual, in the AdS / CFT correspondence large (small) distances in the five- dimensional 
geometry, corresponds to small (large) distances in the gauge theory side The 
variable U represents a distance in the AdS side and it is interpreted as an energy variable 
in the gauge theory side. It is instructive to find the various limits of the metric ( p^ ) and 
of the quark-antiquark potential, by examining the behaviour of the harmonic function 
H as we go from the UV to the IR. The UV corresponds to energies large enough so that 
the vev's can be neglected. The deep IR corresponds to energies very close to the vev's. 

4.1.1 Behaviour for U ^ tq 

In the UV for U ^ vq since the harmonic takes the form 

~ — , for t/ > ro , (75) 
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the metric becomes that for AdS^ x with each factor having radius R in string units, 
Then the original 5*0(6) global symmetry of the metric is restored. That corresponds to 
the fact that in the deep UV the values of the vev's can be neglected and we should be 
describing an SU{Nk) SYM theory. In this regime the integrals can be evaluated and 



give for the potential energy between a quark and antiquark the results of |]T9|. Hence 
we have 

(76) 



This behaviour is similar to the one found from pure gauge theory considerations in (15) 
for the interaction of charges at the same site j = k. 



4.1.2 Behaviour for U — vq ^ ro/N 

For distances from the location of the centers much larger than the distance between the 
different centers we may approximate the harmonic function by 

H = R^ 2 ^ , for f/ - ro > ro/N , (77) 

((^72 + ^2)2 _ 4^2^2)3/2 

since then E/^r ~ 1. Then we see that the discrete part Z^- of the symmetry group 
50(4) X Zn becomes a continuous U{1). Then the resulting integrals for the length and 
the energy are quite complicated (see eqs. (A. 6) and (A. 7) of 0]) and it does not seem 
possible to compute them explicitly, let alone to explicitly obtain the energy as a function 
of the separation length. This approximation is valid for L ^ R\N/rQ. The situation 
here is similar to the one described from a gauge theory point of view by eq. ( P4D (for 
J = k). 

As a separate remark we note that the metric ( p7| ) with the above harmonic (|7^ has 
a ring singularity defined by the equations U"^ = v? = x\ + x1 = r^. It is obvious that 
this is an artifact of the continuous approximation which is not valid near the singularity, 
where the harmonic which should be used is given by ([8^) below and the metric then is 
completely non-singular. The same resolution of singularities applies to solutions that 
have been used for studies of the Coulomb branch of the A/" = 4 SYM gauge theory in, for 
instance, 0, ^ These solutions, take the form (|67D, but unlike (|69D, the distribution 



of D3-branes is continuous instead of discrete. It is obvious that the solutions will become 
non-singular if we replace the continuous by a discrete distribution of D3-branes with the 
given continuous limit. An example of that was already given in the appendix of 
representing the discrete version of a continuous distribution of D3-branes on a disc. 

As another side remark we note that the solution with harmonic ( |69D as well as 
its continuous approximation (|77|), do not have a five- dimensional gauged supergravity 
origin. 
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4.1.3 Behaviour for tq/N <^U — tq <^ tq 



For distances from tlie centers much larger than the distance between the different centers, 
but also much smaller than the radius of the circle, we may approximated the harmonic 
function by 

H ~ , for ro/AT < f/ - tq < tq . (78) 

In this case we have a smeared D3-brane solution along one transverse direction which 
is T-dual to a D4-brane solution. Also the global symmetry has been enhanced from an 
5*0(4) X f/(l) into an 5*0 (5). Since a D4-brane solution corresponds to a five-dimensional 
theory we expect that there will be a behaviour for the quark-antiquark potential. 
Indeed we obtain the result [0 



(79) 



The behaviour is similar to (|25|) from the usual gauge theory point of view. 

4.1.4 Behaviour for U — tq <^ tq/N 

In the deep IR for energy scales very close to the vev's 

H-^, for f/ - ro < ro/N , (80) 
where Rf = Airg^k is the 't Hooft coupling for an SU {k) gauge theory. Then we have 



This is a Coulombic behaviour for a gauge group SU{k) and is similar to that in (|T3D 
from the usual gauge theory point of view. 

4.1.5 Discrete infinite cirray 

Finally, we may also find a regime where the original circular discrete array is approxi- 
mated by a discrete infinite array on a straight line. This is possible for distances from 
the centers much smaller that the radius of the circle. From the gauge theory point of 
view at weak coupling the analogous potential is described by (|26|). The solution has a 
metric given by ( |57[ ) with harmonic 

oo j-,4 

(^2 + _ 2T:vjyf 

Rj f vsmh{U/v) cosh(f//f) cos(x6/t;) - 1 \ , . 



4^3y2 ycosh(f//f ) — cos(a;6/w) (cosh(f//t>) — cos(x6/f))^ 
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where U"^ = x\ + x\ + . . . + x\ and v = vq/N. Translated into separation distances for the 
quark- ant iquark pair, this approximation is vahd when L ^ R\\/N /tq. The trajectory 
that we will use to compute the potential amounts to setting xg = 0. Then we utilize the 
general formulae ( [73| ) and ( [7^ ) with R replaced by i?4 and the function / being given by 

S(V)^..U^[^MUIi2.))^JIM_y, (83, 

As we move from U ^ v to U '^v this case interpolates between those described before 
by ([79| ) and (pT]). It is not possible to compute the energy as a function of the separation 
explicitly, but we may use the above expression in order to find the first correction to the 
behaviour in ([79|). Expanding ( [83| ) foiU^v and keeping the two leading terms 
we find that 



where the Ci and C2 are two constants of oder l.|^ Since 1 <^ ^^5^ -C N the corrections 
to the leading behaviour are small. However, they do not obey a power-like law as 
the corrections to (noted after ([T^)). Similarly, the corrections to (0) are found by 
expanding (|83|) for U <^ v. It is easy to see that these corrections are organized in powers 



of -^j- <^ 1, in contrast to the corrections to ([13|) which are of the Yukawa type (noted 



again below ([lOD ). Hence, it seems that the mechanisms for developing an effective fifth 
dimension are different at weak and strong coupling. 

In will be interesting to further extend the results of this paper in order to incorporate 
fermions and to investigate the renormalization of the gauge coupling constant. Also an 
extension to non-commutative gauge theories could be of some interest. 
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A Varying gauge coupling 

In this appendix we first study differential operators defined on the lattice of the type 

Tj = q]qj, J = 0,1,...,N, (A.l) 

for some operator qj and which have energy eigenvalues that are positive semi-definite, 
i.e. > 0. In the case of the operator (^ discussed in the main text we have that 

^Precisely, ci = Vs^ - l) [tuMT - ^'^^^ ^'^^^ ^2 = 4^ (fiw)' - O-'^^^' 
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1j ~ '^fji^'^^ " 1)5 already noted. If a state 0j_o satisfies the equation q'j</>j,o = 0, then 
automatically we have a solution of Tj^j = corresponding to a zero mass eigenvalue. 
Define now the lattice differential operator Tj obtained by reversing the order of qj and 
q] in 

f,=q,q], j = 0,l,...,N. (A.2) 

The discussion next is completely parallel to that in supersymmetric quantum mechanics 
(SQM) [0. The spectra and eigenfunctions of and of Tj and Tj 

Ml = Ml^,, M^ = 0, 

h,n = J^qj<j)j,n+l , (A.3) 

^hn+l = -^(lj<Pj,n ■ 
Mn 

Hence, the spectra of the operators Tj and Tj are identical except for the zero mode. As 
we have seen, boundary conditions may project it out of the spectrum in which case the 
spectra are identical. Moreover, acting with qj (gj) we may convert an eigenvector of Tj 
{Tj) into an eigenvector of Tj {Tj). 

Let us define the matrices 



--(o'l)' 



It is trivial that these obey the closed superalgbebra 

[H,, Q,] = [H,, Q]] = {g„ Q,} = {Q], Q]} = , {Q,, Q]} = H, . (A.5) 

In terms of this superalgebra, the fact that Qj and commute with Hj is responsible 
for the degeneracy of the spectra. 

As an application of the above general results, let us consider again the low energy 
effective action in the lattice (0), but now with varying gauge coupling gf j and Higgs 
values Vj . The relevant term derived from (|ID corresponding to the gauge field fluctuations 
is 

„ N N-l 

S = -TTly-TT / E t^^A)' + E ^X94,j-,iAi,^' - 94,jAif ■ (A.6) 
Consider first the case of constant vev's Vj = v. The mass term above can be written as 

N 

i=o 

where the operator Tj is defined as 

Tj = v'g^jie-''^ - l)(e'^^ - 1)^74,, = -"^v^Qaj smh\d,/2)g^j . (A.8) 
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Hence, the operator Tj has precisely the form ( [A.2| ) with qj = gij{e~'^^ ^1)- Using the 
formahsm that we developed we may map this problem into a problem for the super- 
partner operator Tj in ( |A.lj ). Comparing with the case of constant gauge coupling and 
varying vev's considered in section 2, we find that they are identical after we identify 
v'^q] . , 1 by 1 . 

As a specific application consider the case of couplings g^j = gi\/], j > 1- After 
using the results of subsection 3.3 we find that the set of eigenf unctions of the operator 
Tj are 

0,- ~ LA IV?) , m = — . (A. 9) 

94V 



Using ( |A.3| ) we find that the eigenfunctions of the operator ( |A.8| ) are 



j(L,_i(m2)-L,(m2)) . (A.IO) 

We may easily verify, using recursion relations for the Laguerre polynomials, that this 
solution satisfies the eigenvalue equation for the operator Tj in ( |A.8[ ). We have left un- 
determined the overall normalization constant, since this would depend on the boundary 
conditions that we will impose. Note that a choice of a boundary condition for either (pj 
or (pj implies the boundary condition for the other, if we want that the mass spectra are 
also related according to ( |A.3|) . Of course we may still impose independent boundary 



conditions in which case the two spectra are not related. 

Finally, we note that in the most general case of varying vev's and gauge couplings, 
the operator Tj still takes the form ( [A.l|) with qj = Vj{e'^^ — l)5'4,j. 
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